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Abstract 

Perturbative subcritical series expansion for the steady properties of a class of one-dimensional 
nonequilibrium models characterized by multiple-reaction rules are presented here. We developed 
long series expansions for three nonequilibrium models: the pair-creation contact process, the 
A-pair-creation contact process, which is closely related system to the previous model, and the 
triplet-creation contact process. The long series allowed us to obtain accurate estimates for the 
critical point and critical exponents. Numerical simulations are also performed and compared with 
the series expansions results. 

PACS numbers: 02.50.Ga,05.10.-a,05.70.Ln 
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I. INTRODUCTION 



Nonequilibrium systems have been used to describe a variety of problems in physics, 
chemistry, biology and other areas. Different classes of nonequilibrium systems have been 
studied and a special attention has been devoted to systems with absorbing states. The 
most studied system with absorbing states is the contact process (CP) The basic CP is 
composed of spontaneous annihilation of particles and creation of particles in empty sites 
provided they have at least one nearest neighbor site occupied by a particle. The criti- 
cal behavior of CP belongs to the directed percolation (DP) universality class P|. Several 
approaches have been aplied for describing the CP, such as numerical simulations 2 L con- 
tinuous description by means of a Langevin equation 0, ^, renormalization group js, y| and 
series expansions 

mm. 

Very accurate estimates for the critical point and critical 
exponents have been obtained for the basic CP model. However for the models that concern 
us here, the avaiable results come from only numerical simulations. The use of different 
tecniques may be a useful tool for studying the behavior of systems whose avaiable results 
are controversial. 

In this paper, we developed subcritical perturbative series expansions for the steady state 
properties of a class of one- dimensional nonequilibrium models characterized by catalytic 
creation of particles in the presence of m-mers. We have considered here creation of particles 
in the presence of pairs of particles (m = 2) and triplet of particles (m = 3). To exemplify, we 
developed the series expansion for three nonequilibrium models: the pair-creation contact 
process (PCCP) jlll . 11^ . the A-pair-creation contact process (APCCP), which is closely 
related to the PCCP, except that an empty site in the presence of at least a pair of particles 
becomes occupied with rate A (in the PCCP it becomes occupied with rate A/2 times the 
number of pairs of adjacent particles), and the triplet-creation contact process (TCCP)^, 



12| |. Very precise estimates of the critical behavior and critical exponents are obtained 



from the analysis of the series. We compared our results with their respective numerical 
simulations. 
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II. OPERATOR FORMALISM 



Let us consider a one-dimensional lattice with N sites. The system envolves in the time 
according to a Markovian process with local and irreversible rules. The time evolution of 
the probability P{r],t) of a given configuration rj = (?7i, ?72, ?7iv) is given by the master 
equation, 

d ^ 

-P{7j, t) = J2{Mv')Piv\ t) - W,{7j)P{7J, t)}, (1) 

i 

where 77' = (771,772, 1 — Vii ■■■iVn)- The total rate Wi{r]) of the interacting models studied 
here is composed of two parts 

w,(77) =<(r7)+<(r/), (2) 

where the rates wf{r]) and (77) take into account the annihilation and the catalytic creation 
of particles, respectively. The annihilation of a single particle is illustrated by the scheme 
• — o, and the catalytic creation of particles hj o • • ••• and o • •• • • ••, for 

the PCCP and TCCP, respectively. The scheme for the APCCP is similar to the PCCP, 
but here if an empty site has at least one pair of adjacent particles, a new particle will be 
created with rate A. More precisely these rules are given by 

<{V)-Vi: (3) 

for the annihilation subprocess, and w^(?7) by 

^iiv) = ~ Vi){Vi+lVi+2 + Vi-2Vi-l): (4) 

for the PCCP, 

tJ^iiv) = ^(1 - ?7j)(^i+i^j+2 + ?7i-2?7i-i - r]i+i'ni+2r]i-2'ni-i), (5) 
for the APCCP, and 

^ ~ Vi){Vi+lVi+2Vi+3 + 77i-3?7i-2?7i-l), (6) 

for the TCCP. 

Before developing the series expansion, it is necessary to write down the master equation 
in terms of creation and annihilation operators. The base states corresponding to a given 
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site i of the lattice are {rji) with \rii) = |o) or \rii) = |«) according to whether site i is vacant 
or occupied by a particle, respectively. The creation and annihilation operators for the site 
i are defined in the following manner 

A+|r/,) = (l-r7,)|l-r/,), (7) 

and 

A\Vi) = Vi\'^-Vi): (8) 

and they satisfy the property AfAi + AiAf = 1. 

Introducing the probability vector defined by 

\^{t))=j:piv,m, (9) 

V 

where \r]) = Ylf=i0\r]i) = \r]i,r]2, ...,r]j\f) is the vector defined by the direct product of the 
base vectors. Substituting Eq. (0) in to Eq. (P) and using Eqs. (|7j) and (jHI), the time 
evolution for the probability vector is given by, 

j^m)) = w\m), (10) 

where the operator W = Wq + XV is a sum of the unperturbed term Wo and a perturbed 
term AV^. The operator Wo, that takes into account only the annihilation subprocess, is a 
nontinteraction term, given by 

Wo = Y.iA-^tA). (11) 

i 

Each term of the summation has the following set of right and left eigenvectors 

|0) = |o), (0| = (o| + (.|, (12) 

with eigenvalue Aq = and 

|l)^-|o) + |.), (13) 

with eigenvalue Ai = —1. The operator V, correspondint to the catalytic creation of parti- 
cles, is an interacting term, given by 

V = ^Y1 Si{ni+ini+2 + ni^2ni-i), (14) 

i 



for the PCCP, 

V = Y1 Si{ni+ini+2 + ni_2ni-i - ni+ini+2ni-2ni-i) , (15) 

i 

for the APCCP, and 

^ = ^ XI 5'j(nj+ini+2ni+3 + ni_3ni_2?2i_i), (16) 

i 

for the TCCP, where the operator Si is given by Si = Af — AiAf and rii = AfAi is the 
operator number. 

To find the steady vector lip), that satisfies the steady condition (Wq + \V)\ip) = 0, we 
assume that 

oo 

l^) = lV'o) + E^'lV'.), (17) 

where \ipo) is the steady solution of the non-interacting term Wo satisfying the stationary 
condition 

Woltpo) = 0. (18) 

T 



a 



le vectors ^pi) can be generated recursively from the initial state |V'o)- Following Dickman 
7 1, we get the following recursion relation 

1^,) = -RVliJe-i). (19) 



The operator R is the inverse of Wq in the subspace of vectors with nonzero eigenvalues and 
given by 

^ = E \<Pn)^{<Pnl (20) 
n(^0) 

where and are right and left eigenvectors of Wq, respectively, with nonzero eigen- 
value An- 

We notice that the steady solution of the noninteracting operator Wq corresponds to the 
vacuum \iIjq) = |0o) = |-0.). Since the creation of particles is catalytic, then if we start from 
the vacuum state, we will obtain a trivial steady vector namely = \ipo). To overcome 
this problem, it is necessary to introduce a modification on the rules of the models. The 
necessity of introducing a small modification on systems with absorging states in order to 
get nontrivial steady states has been considered previously by Tome and de Oliveira 



and by de Oliveira 



li 
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III. GENERATING THE SUBCRITICAL SERIES 



The modification we have made consists in introducing a spontaneous creation of particles 
in two specified adjacent sites for the PCCP and APCCP. The chosen sites are i = and 
i = 1, so that the rates Wq^t]) and wlirf) are changed to 

<(r/) = (l-g)r/o + g(l-r7o), (21) 

and 

<(r7) = (l-g)r/i + g(l-r7i), (22) 

where q is supposed to be a small parameter. This modification leads to the following 
expression to the operator Wq 

Wo = Y. W^o^ + g(5oo + 5oi - Wo, - Woi). (23) 

i 

The steady state of Wq is not the vacuum state anymore. Now, it is given by 

|^/.o) = |.0.) + 2g|.10.) + g2|.ii.), (24) 

where all sites before and after the symbol "." are empty. 

Two remarks are in order. First, only the last term in l^/^o) will give nonzero contributions 
to the expansion so that l^/'^), i > 1, will be of the order g^. Second, although the change 
in Wq will cause a change in R, only the terms of zero order in the expansion in q, given by 
the right-hand side of Eq. ((201), will necessary since the corrections in R will contribute 
to terms of order larger than g^. For instance, the two first vectors, {ipi) and IV'2); for the 
PCCP are given by 

= g2{2|.l.) + + 1.101.) + ^1.111.)}, (25) 

and 

l^2) = g'{^|.i.) + ^|.ii.) + ^|.ioi.) + ^|.iii.) + ^|.iooi.) + ^|.iioi.) + ^|^ 

(26) 

The translational invariance of the system is assumed. 

For the TCCP, the rates w'^{ri), i = 0,1,2 are modified similarly and an analogous initial 
vector \ipo) is obtained. However, the vectors lipi), i > I, will be of the order q^. 
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The series expansions for the probabihty vector \ip) obtained here are equivalent to the 
Laplace transform of the time dependent vector probability \^{t)) in the subcritical 

regime. If we assume that |\E'(s)) can be expanded in powers of A, 

\^{s)) = l^o) + A|^i) + A^l^s) + • • • , (27) 

where 

\'^o) = {s-Wor'\Xo), (28) 

and 

= {s-Wo)-'V\^i-i), (29) 

where |Xo) = | • •) for the PCCP and APCCP and |Xo) = !•••) for the TCCP. The two 
first vectors for the PCCP is given by 

l^o) = ^1-00.) + 2si|.10.) + s2l.il.), (30) 

and 

1^1) = 2s2(si|.l.) + s2l.ll.) + S2I.IOI.) + S3I.III.)), (31) 

where s,. = l/(s + r). In the limit s ^ 0, Eq. (jHT|) becomes identical (by a factor 2g^) to 
the Eq. (f^3j) . that is |\E'i) = |'?/'i)/2g^. The next orders of the expansion will also produce 
vectors that follows a similar relationship, namely |\E'£) = \ip£)/2q'^. Therefore, the steady- 
state vector |\E') has a close relationship with the Laplace transform |^E'(s)) of the time 
dependent vector probability {"^{t)) in the subcritical regime. 

IV. ANALYSIS OF THE SERIES 

To calculate the coefficients of |\&£) in the base \7]) we have built a computational algorithm 
to take account of all configurations. The configuration can be expressed in terms of a binary 
number rji + 7^22 + 7732^ + . . . representing the vector \ri). For example, the binary number 
1101 corresponds to the configuration |.1101.) and we need to store only the value of the 
coefficient of 1101. By this procedure we were able to determine the coefficients of all vectors 
I'^i) up to the 26th order in A for the PCCP (and APCCP) and to the 25th order for the 
TCCP. 

From the series expansion of the vector lip), it is possible to determine several quantities, 
such as survival probability, the total number of particles, and the correlation function. In 
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TABLE I: Coefficients for tlie series expansion for total number of particles A^" corresponding to 
tfie PCCP, APCCP, and TCCP 
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TABLE II: Continued from Tabled 
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this paper, however, we will be concerned only with the series expansion for the total number 
of particles A^, given by 

N={.0.\J2n,\^p). (32) 

i 

One can show that the coefficient of in the expansion for is simply the coefficient of 
|.l.) in \ip£). This allows us to get a longer series for the number of particles. For the PCCP 
and APCCP we obtained 38 terms and for the TCCP we obtained 33 terms. The resulting 
series for the total number of particles of the three models considered here are listed in Table 

m 

From the series expansion of a given quantity, in the present case, A^, we can determine 
the critical point and its corresponding critical exponent by means of a Pade analysis. Since 
the series developed here is related to the Laplace transform of the total number of particles. 



both will have the same critical behavior namely 



AT - (A- Ae)-"ii(i+''), 



(33) 



where and rj are the exponents related to the time correlation length and to the growth 
of the number of particles, respectively. 

A preliminar analysis is done by performing unbiased estimates for determining both the 
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FIG. 1: Biased estimates of = z^||(l + as a function of Ac derived from the Fade approximants 
to the series (A — Xc){d/dX) InN = 6 evaluated at Ac for the FCCF. The approximants shown are 
[17/18], [18/17], [18,19], [19/18], and [18/18]. 



critical point Ac and the critical exponent by means of the Fade aproximants [l^lla]. This 
approach consists of analysing the serie (d/dX) InN by a Fade approximant. The critical 
exponent and the critical parameter Ac are obtained from the pole and the residue at this 
pole, respectively. We have obtained unbiased analysis for the three models considered 
here. However, they give us estimates that does not seem to improve significatively when 
we consider higher-order Fade approximants. For example, for the FCCF the approximant 
[13/13] gives Ac = 7.62 and z/||(l+?7) = 2.71 whereas the approximant [16/16] gives Ac = 7.54 
and z/||(l + rj) = 2.56. 



Much more reliable estimates are obtained when we perform biased analysis, w. 



lich is set 



ISlQ- For 



up by looking at Fade aproximants to the series (A — Xc){d/dX) InP = 6 j l7l . 
a trial value of Ac, we develop the serie above obtaining ^(Ac) for a given Fade aproximant 
[m/n]. We can build curves for different Fade approximants by repeating this procedure for 
several trials Ac and we expect that they intercept at the critical point (Ac, 6'(Ac)). In the 
Figs. HJ El and El we plotted the curves obtained by considering different Fade aproximants 
for the three models. 

From the Figs. Q El and we see a very narrow intersection of the Fade approximants. 
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FIG. 2: Biased estimates of ^ = z^||(l + as a function of Ac derived from the Fade approximants 
to the series (A — Xc){d/dX) InA^ = 9 evaluated at Ac for the AFCCF. The approximants shown 
are [18/18], [17,18], [18/19], and [19/18]. 



TABLE III: Biased estimates for Ac and z^||(l + 7y) from the Fade approximants for the three models 
considered here together with the values for the basic contact process (CF) 



Model 


Ac 


i^||(l + r/) 


FCCF 


7.4650(6) 


2.274(3) 


AFCCF 


3.9553(5) 


2.272(4) 


TCCF 


12.01(2) 


2.26(2) 


CF 


3.29782 


2.2772 



revealing the utility of this approach. However, as pointed out by Guttmann it is 
difficult to estimate uncertainties in series calculations. Thus, in order to give a more 
realistic estimate of the quantities measured here and their associated uncertainties, we 
have estimated them by taking into account the first and last crossings among various Pade 
approximants. The values of the critical parameters obtained for the three models are 
summarized in the Table HVl The estimates of Ac for the PCCP and TCCP are in excellent 
agreement with the corres pon ding values Ac = 7.464(2) and Ac = 12.00(1) obtained from 
numerical simulations 



jrrespon 

HQ. 
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FIG. 3: Biased estimates oi 9 = v ^^{1 + rj) as a function of Ac derived from the Fade approximants 
to the series (A — \c){d/dX) InN = 9 evaluated at Ac for the TCCF. The approximants shown are 
[15/15], [15,16], [16,15], and [16,16]. 




1 I — ' — ' — ' — ' — ' — ' — ' — ' — ' — I 

2 4 6 8 10 

t 

FIG. 4: Log-log plot of the number of particles Ng versus the time t for some values of A for the 
APCCP. From top to bottom A = 3.97, 3.9558, 3.9553, 3.95, and 3.93. 

V. NUMERICAL SIMULATION 

As a check of the accuracy of the results obtained here, we have performed spreading 
simulations for the APCCP, since its critical point is unknown in the literature. Following 
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Grassberger and de la Torre [20|, we started from a initial configuration close to the ab- 
sorbing state with only two adjacent particles, we can study the time evolution the survival 
probablility Ps{t), the mean number of particles Ns{t), and the mean-square distance of 
the particles from the origin R{t). At the critical point, these quantities are governed by 
power-laws whose their related critical exponents are named S, rj and z, respectively. Off- 
critical point, we expect deviations from the power-law behavior. In the Fig. |3J we plotted 
the quantity Ns{t) versus the time t for some values of A. Analogous analysis can be done 
for determing the exponents 6 and z. At the critical value Ac = 3.9553, our data for the 
three quantities Ns(t), Ps(t), and R(t) follow indeed a power law behavior, whose critical 
exponents are consistent with those of the DP universality class. 



VI. CONCLUSION 



We have derived subcritical series expansions for studying the critical behavior of three 
nonequilibrium systems characterized by multiple-creation of particles. Although series ex- 
pansion have been applied sucessfully for the contact process and similar models [sl lol Il7|. 
this is the first time that these systems, with multi-reaction rules, has been treated by means 
of a technique other than numerical simulations. With exception of the TCCP, whose value 
of Ac are in the same level of precision of numerical simulation estimates, the subcritical 
series expansion give us the best estimates for the critical point of the models considered 
here. The critical exponents are consistent with those related to models beloging in the 
DP universality class. We remark finally that the present approach may be very useful to 
determine the critical behavior and universality classes for other nonequilibrium systems. 
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